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_C ■ Abstract. We apply the method of skew-orthogonal polynomials (SOP) in the 

C^ , complex plane to asymmetric random matrices with real elements, belonging to two 

different classes. Explicit integral representations valid for arbitrary weight functions 

are derived for the SOP and for their Cauchy transforms, given as expectation values 

of traces and determinants or their inverses, respectively. Our proof uses the fact 

^ ' that the joint probability distribution function for all combinations of real eigenvalues 

C*) , and complex conjugate eigenvalue pairs can be written as a product. Examples for the 

SOP arc given in terms of Laguerre polynomials for the chiral ensemble (also called the 



£s^i ■ non-Hcrmitian real Wishart-Laguerre ensemble), both without and with the insertion 

of characteristic polynomials. Such characteristic polynomials play the role of mass 
f^i ' terms in applications to complex Dirac spectra in field theory. In addition, for the 

^D . elliptic real Ginibre ensemble we recover the SOP of Forrester and Nagao in terms of 

Hermite polynomials. 
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1. Introduction 

Classical orthogonal polynomials (OP) are one of the principal standard tools used to 
solve problems in Random Matrix Theory (RMT). The three classical Wigner-Dyson 
ensembles with Gaussian elements can be solved in terms of Hermite polynomials, 
whereas their chiral counterparts require the use of Laguerre polynomials. Whilst for 
the symmetry classes with unitary invariance (/3 = 2) the corresponding scalar product 
is symmetric with the standard weight, for the classes with orthogonal (/3 = 1) or 
symplectic ((3 = 4) symmetry the scalar product becomes skew-symmetric, with the 
details - including the weight - dependent on the symmetry class. The corresponding 
polynomials are then called skew-orthogonal polynomials (SOP). For details of all these 
cases we refer to [1], as well as to [2] and [3] for reviews on SOP. What all these ensembles 
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have in common is that each of their solutions can be expressed using the kernel of the 
corresponding (S)OP as a building block. 

The solution of RMT in terms of (S)OP is exact for finite matrix size N. Moreover, 
when taking one of the possible large- iV limits in the bulk of the spectrum, or at the 
soft or hard edge, the standard Plancherel-Rotach asymptotics can be used. Much work 
has been done on the question of universality, i.e. the extent to which the asymptotics 
also hold for non-Gaussian weights, see e.g. [I] for a review. 

This setup of (S)OP has been generalised to the complex plane in order to solve 
non-Hermitian RMT. However, for the standard Ginibre ensembles the (S)OP are 
simply given in terms of monic powers (which holds for all weights that are rotationally 
invariant in C). Only when considering so-called elliptic deformations of the Gaussian 
Ginibre weight do Hermite polynomials on C appear, as was first observed in [5]. The 
corresponding quaternionic (/3 = 4) elliptic Ginibre ensemble was solved in terms of some 
Hermite SOP in [6], and only very recently was the real elliptic Ginibre case (/3 = 1) 
solved in terms of another set of Hermite SOP [7j. 

The chiral counterparts of two of these ensembles were introduced in [8] (/3 = 2) and 
[9] (P = 4), where they were solved in terms of Laguerre OP on C and Laguerre SOP 
respectively. The obvious question then arises about the existence of Laguerre SOP for 
the chiral /3 = 1 ensemble [TUJ [II] which we answer affirmatively in this article, thereby 
completing the set of classical Hermite and Laguerre (S)OP in the complex plane for 
these non-Hermitian RMT. 

One complication arises for the ensembles with Laguerre (S)OP in the complex 
plane: due to the integration over angular variables the elliptic ensembles that are 
Gaussian in terms of the matrix elements lead to non-Gaussian weight functions for the 
eigenvalues (for (3 = 1 elliptic Ginibre we have also a complementary error function). 
The Bessel function of the second kind appearing here for all three /3 = 1,2,4 makes 
the orthogonality question much more involved. 

For this reason we first provide a new integral representation for the = 1 SOP, 
valid in both the elliptic Ginibre and chiral symmetry classes for arbitrary weight 
functions. In view of earlier results for SOP on R for both j3 = 1, 4 [i~2| IT3] . as well as 
for complex SOP for /3 = 4 (6j [9], this representation comes very naturally. Moreover, 
it was shown very recently in [2] that the /3 = 1 and /3 = 4 Ginibre ensembles can 
be treated on an equal footing. We will rederive this relation amongst these symmetry 
classes from a different angle. We will also derive a new integral representation for the 
Cauchy transforms of the SOP on C valid for both /3 = 1 and 4. This extends the 
expression for Cauchy transforms on C for (3 = 2 in |15j . 

One important ingredient necessary in order to derive these results is the 
factorisation of the joint probability distribution function (jpdf) for /3 = 1, which is 
originally given by a sum over all possible combinations of real and complex conjugate 
eigenvalues [T6| [T7] . Such a factorisation, which uses the symmetrisation over all 
eigenvalues, might have been expected from the fact that the partition function can 
be written as a single Pfaffian over double integrals pi 
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Our integral representation allows us to derive the /3 = 1 Laguerre SOP on C in a 
straightforward fashion. The known Hermite SOP on C of [7] also follow easily. As a 
third and important example for our general formalism we explicitly compute the SOP 
for the chiral /3 = 1 ensemble with the insertion of mass terms. Such insertions play 
a crucial role in the application of RMT to the complex Dirac operator spectrum in 
Quantum Chromodynamics (QCD) and related field theories at small quark chemical 
potential in the low density phase, see e.g. [19] for reviews, as well as very recently in 
the high density phase corresponding to maximal non-Hermiticity |20j . 

For this third example we exploit the fact that averages (and ratios) of the required 
characteristic polynomials have been computed very recently for the non-Hermitian 
/3 — 1 symmetry classes in |21j . Together with our interpretation of the building blocks 
there as the kernel, SOP and their Cauchy transforms this completes the analogy to 
earlier computations of such averages and ratios for /3 = 2 in [221 E3 E3] and /3 = 4 in 

js |23]. 

This article is organised as follows. In section |2] we recall the definition of the 
two RMT with real asymmetric matrices including mass terms, their respective weight 
functions and eigenvalue representations. Section [3] is devoted to a factorisation proof 
of the jpdf of real and complex eigenvalues, where we give two different arguments. The 
new integral representations for the SOP and their Cauchy transforms are then shown 
in subsections 14.11 and 14.21 respectively. In section |5] we provide three explicit examples 
for SOP including Hermite in subsection 15.11 Laguerre SOP in subsection 15.21 and SOP 
including masses in subsection 15.31 In the Appendices we collect together short proofs 
for some mathematical identities used in the text. 

2. The Matrix Models 

We will show how to solve the following two matrix models of real asymmetric matrices 
in terms of skew-orthogonal polynomials (SOP) in the complex plane. 

The first model is given by the chiral extension of the elliptic real Ginibre ensemble 



A' 



/ 



*£"'({".» ~ /_ dAdBj[det[ ""•£*" ,„,„,„ ) (2-1) 

1 



R 2iY(JV+„) " y D l'lfn-(N+v)y.(N+i>) 



x exp 



--?7 + Tr(A4 T + BB T ) + ?7_Tr(AB 3 



with r)± = —^- . (2.2) 



l±/i 2 
4/x 2 

Here A and B are real asymmetric matrices of size N x (N + u), and /i 6 (0, 1] is 
a non-Hermiticity parameter. The integration runs over all independent real matrix 
elements of A and B with a flat measure. The product of determinants or characteristic 
polynomials is motivated by the addition of Nf quark flavours in applications to QCD at 
finite density [191 I2Q]- The model can be written as a Gaussian two-matrix model with 
two independent real asymmetric matrices P and Q, with A = P + fiQ and B = P — \xQ 
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(see PH]). In the limit \x — > the model reduces to the chiral Gaussian Orthogonal 
Ensemble (chGOE). 

The second matrix model is also a generalisation of the elliptic real Ginibre ensemble 
and is given by 

N f 



1-T 



Tr( J J 1 - tJ 2 ) 



(2.3) 



42\{™})~ I a dJn|det[J-<m,]| a exp 

Here J is a real asymmetric matrix of size N 2 , and r £ [0, 1) is the non-Hermiticity 
parameter. We again integrate over all the independent matrix elements of J. The model 
can alternatively be written as a two-matrix model with symmetric and anti-symmetric 
matrices S and A with Gaussian elements, where J = S + A\/(l — t)/(1 + r), see also 
[25] for Nf = 0. In the limit r — > 1 the model reduces to the GOE when Nf = 0. The 
extra determinants correspond to the imaginary mass terms coming in pairs of opposite 
sign. These are also motivated from applications to QCD with a chemical potential, but 
this time in three dimensions, see [26] . 

Once we switch to an eigenvalue basis for the Dirac matrix T> for the first model, 
where 



<W A 



B T Or 



and to the eigenvalues of the matrix J for the second model, both models can be treated 
along the same lines. Because the characteristic equation for both T> and J is real its 
solutions are either real or come in complex conjugate pairs. However, because of the 
chirality of T> there is a peculiarity here: 

N 

= det[Al 27V+ , -V] = A v det[A 2 ljv - AB T ] = A" J^[(A 2 - A 2 ) . (2.5) 

i=i 

Whilst the A 2 are indeed either real or come in complex conjugate pairs, the Dirac 
eigenvalues Aj are consequently real (A| > 0), purely imaginary (A 2 < 0), or come in 
quadruplets (±Aj,±A*); there are also v generic zero-eigenvalues. For simplicity and 
to keep the presentation of the two models parallel we will mainly consider changed 
variables Zj = A 2 in the following. 

Following [271 EH] the partition functions in eqs. (12.11) and (12.31) above can be written 
as follows, where the normalisation is to be determined later, see eq. ( 14.161) . 

N 2n+x r N-n „ 

Z2N+ x = N] -^2 II / dXk II / d 2 ZmP2n+ X ,N-n(x,Z,Z*) . (2.6) 

n=0 fc=l jj m=l „ 

Here we sum over all the possible ways of splitting the total number (2A^ + x) of 
eigenvalues into K = In + x rea I eigenvalues {xk} and M = N — n complex conjugate 
eigenvalue pairs {z m , z^}. A product with an upper limit less than its lower limit is 
defined as unity. Note that in this expression we have only one complex integration for 
each complex conjugate eigenvalue pair. The differentials of the complex eigenvalues are 
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defined over the real and imaginary part, i.e. d 2 z m = d$lez m d^$mz m . In the following 
we treat the cases of an even (x = 0) and odd (x = 1) total number of eigenvalues on 
the same footing. The jpdf for a fixed number K of real eigenvalues and M complex 
eigenvalue pairs is defined as 

K M 

Pk,m(x,z,z*) = ] \h{x k ) ] [ [g(z m ,z* m )2iQ^sm.z m )\ A K+2M ({#}, {z, z*}) 



fc=l m=l 

K 

X 

fc=2 m=2 



H<d(x k - a*_i) J] Q(^ez m - Uez m ^) (2.7) 



with the weight specified in eqs. (J2.9P and (12.111) below [j. Here is the Heaviside 
distribution, and the Vandermonde determinant is defined as 



JV 

A*({*}) = II(^ - ^ = , /et [4- 1 ] . (2.8) 

-*--*- Ka,b<N L 

In eq. ( 12. 7ft we explicitly specify the number X of real eigenvalues and M complex 
eigenvalue pairs, with the set of arguments labelled as Xi, . . . , x^, zi, z*, . . . , zm, z*m m 
Ak+2m({x}, { z ) z *})- The factors 2i®{^svaz rn ) and the ordering of the real eigenvalues 
0(xfc — Xk-i) in eq. (12.71) allow us to omit the modulus sign around the Vandermonde 
determinant. The ordering of the real parts G(3?e z m — 3f?e2 m _i) is needed to make the 
transformation to upper triangular 2x2 block form of the matrices A and B T (or J) 
unique when computing the Jacobian [11] (see also [28]). The latter and part of the 
former can be dropped later due to the symmetrising integration as will be shown in 
the next section. 

We also mention that the partition function can be written as a single Pfaffian 
[TBI 127] . eq. ( 13. 8ft below, and we come back to the consequences for factorisation in the 
next section. 

We now give the weight functions for our two models. Looking at eq. ( 12. 5ft for the 
chiral model we are only interested in the eigenvalues of the matrix C = AB T . Because 
the Nf extra mass terms compared with [UJ depend only on C their addition to the 
jpdf in [UJ is trivial and so we only give the result. The corresponding weight functions 
in eq. (J2.7P for the real eigenvalues x and complex eigenvalues z = x + iy read 

h ch (x) =2\x\ u l 2 K v/2 {r 1+ \x\)eMV-A\{^ + m)) , (2-9) 

/=i 



gch(zi,z 2 ) = 2\z 1 z 2 \ y/2 eyLv[ri-{z 1 + z 2 )]\\{z 1 + m 2 f )(z 2 + mj) 

/=i 

/°° dt r 1 1 / \ 

— exp - r]\t{z\ + z\) - — K u/2 i^r^tz^ erfc I r)+Vi\z 2 - z^j 



| In contrast to [TT] we distinguish the weights for real and complex eigenvalues by different symbols 
(h and g respectively). 



Skew- OP for chiral real asymmetric random matrices 6 

In addition we have a trivial overall factor arising from the generic zero eigenvalues in 
the mass terms: 

N f 

P% M (x, z, z*) = n m / P k,m{x, z, z*) (2.10) 

/=i 
For the second model eq. (12.31) we have instead 

h Gin (x) = exp[-x 2 ] jj(x 2 + m}) , (2.11) 

/=i 

(\z -z\\ Nf 
gGin(zi, z 2 ) = exp[-z 2 - zl\ erfc I x 1 J \\(z{ + m 2 f )(zj + mj) . 

For the two Gaussian examples above the following relation is satisfied 

lim g{z,z*) = h^tez) 2 , (2.12) 

Smz-tO 

relating the two weights. 

As a general remark here and in the following we can allow for more general weight 
functions h(x) and g(zi,z 2 ) in eq. (12.71) that do not necessarily follow from a matrix 
representation. For example, we could generalise the weights in eqs. (I2.9P and (12. lip by 
multiplying by a factor exp[— V(zi, z 2 )} where V is a polynomial in z\ and z 2 [§|. Moreover 
one can independently choose h and g instead fulfilling the relation (I2.12p . 

3. Factorisation of the joint probability distribution 

In this section we prove that the jpdf inside the partition function can be written in 
a factorised form. For this to be possible, it is essential that we integrate over all the 
eigenvalues, leading to a symmetrisation. Hence this applies equally to the expectation 
value of any operator symmetric in all variables. Examples for this are the computation 
of the gap probability or expectation values leading to integral representations of the 
SOP in the next section. However, such a factorisation can also be found for the /c-point 
density correlation functions when summing over all possibilities of splitting k into real 
and complex eigenvalue pairs. 

Let us first state the result for the partition function eq. ( 12 .6p in terms of a single 
product for the weights 

„ 2N r N 

Z 2N+X = dy x h x (y)Y[ d 2 z k Y[F(z 2j ^ u z 2j ) A x+2N (y,{z}) , (3.1) 

R k=1 C • 7=1 

where we define the anti-symmetric function 

F(z u z 2 ) = ig( Zl , z 2 ) (6(Sm Zl ) - 9(Sm z 2 )) 5 2 (z 2 - z{) (3.2) 

H — h(zi)h(z 2 )SC^m zi)5($sm z 2 )sgn(JRe z 2 — $tez±) . 

§ A so-called harmonic potential could be realised as a matrix model, by multiplying eqs. (|2.ip and 
(|2.3p with cxp[— TrV(AB T )] or cxp[— TrV(J 2 )]. Although at finite N these are formal expressions due 
to lack of convergence this can be dealt with in the large- TV limit. 
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Note that for an even number of variables (% = in eq. ( 13.1ft ) the integration and 
weight for the real variable y have to be dropped, as well as the argument y inside the 
Vandermonde determinant. In eq. ( 13. ip we integrate over 2N independent complex 
variables in contrast to the ordered integration in complex conjugated pairs. The 
standard two-dimensional delta function in eq. ( 13. 2ft reads S 2 (z) = S(x)S(y) for z = x+iy. 

We will prove this factorisation in two different ways. One is by explicitly summing 
up all the terms in eq. (12. 6p to make a single product. The results for expectation 
values of characteristic polynomials in [21] are built up on this idea, although a proof of 
this was not given. The second way is by starting from a single Pfaffian representation 
of the partition function derived in [181 EZ] and using (in reverse) a proof of a slightly 
generalised version of the de Bruijn integral formula. Because both derivations are short 
and illustrate different aspects we decided to present both. 

The first derivation of the factorisation of the jpdf goes as follows. From eqs. (12. 7ft 
and ( 12. 6 P we obtain 

N 2n+ X „ 2n+ X 

z 2N+x = n\ j2 ( N _ n \\ n / dxi h ^ n e ^ ~ x ^ 

n=0 ^ >' 1=1 | j=2 

N-n „ 
X II (d 2z m9{z m ,z^ n )2iQ(Qmz m )) A {2 n+ X )+2(N~n){{x},{z,Z*}) , (3.3) 

m=l c 

where we use that the integrand is totally symmetric under a permutation of two pairs 
of complex conjugated eigenvalues. Dropping the ordering of the real parts leads to a 
factor l/(iV-n)!. 

The ordering of the real variables can be simplified by applying the method of 
integration over alternating variables [1] twice. Pulling the integrations over odd 
variables in, dropping the symmetrisation giving a factor 1/n! and then pulling them 
back out leads to the following result: 

N ^.| 2n+x j. n 

Z 2N+X = Y, nUN-nV 11 / dxi h ^ II & ^ ~ x ^- 1 ) 

n=0 ' K >' 1=1 ^ j=l 

N-n „ 

X II (d 2z mg{z m ,z^ n )2iQ(Qmz rn )) A {2n+x)+2 ( N ^{{x},{z,z*}) . (3.4) 

m=l c 

Isolating the integration dx x we can do the sum over multiple integrations, using the 
binomial formula and the permutation invariance of the integrand under exchanging 
pairs of complex numbers: 



dx\h{x\) I dx 2 h(x 2 ) Q(x2 — x±) 

d 2 z l h{z l ) ( d 2 z 2 h(z 2 )e(^ez 2 -^ez 1 )5(Qmz 1 )5(Qmz 2 ) (3.5) 

c 

b= I d 2 zg(z,z*)2ie(%mz) = J d 2 z x fd 2 z 2 g(z 1 , z 2 ) 2z6(3m^) 5 2 (z 2 - z\) , (3.6) 



c c 
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with 



JV 

E 

n=0 



AH 



niN—n 



n\(N -n)\ 



a n b 



(a + by . 



(3.7) 



Pulling out all 2N independent complex integrations of the iV-fold product this leads 
to eq. (13.11) as claimed, after making the function F(zi, z 2 ) manifestly anti-symmetric. 
We now come to our second argument, starting from the result derived in [T8J [27] . 
This states that, including normalisation, 



Z 



2iV 



N\ Pf 

Ka,b<2N 



I ' d 2 z l d 2 z 2 F{z ll z 2 )[z a l - l z, 



O-li-l n-li-1 



Zn Z-i 



, (3.8) 



where for simplicity we only state the even case, i.e. with 2N eigenvalues. The sign of 

1 



the Pfafnan is defined as in [T] such that the Pfaffian of the matrix 



-1 



<8> Ijv is 



unity. Here F(z\,z 2 ) is the function from eq. ( 13. 2p . By using a slight generalisation of 
the proof of the de Bruijn integral formula, in reverse, 

2N „ N 

[ / d 2 z k T\F(z 2j -i,Z2j) det [{f a (z 2 b-l),ga(Z2b}}] 



k=\ 



i=i 



N\ Pf 

Ka,b<2N 



d 2 ud 2 vF(u,v)[f a (u)g b (v) - f b (u)g a (v)} 



[r- 



(3.9) 



where we refer to Appendix A for a derivation (cf. Appendix C.2 in [29|), we obtain 



2A' 



A 1 



Z 2N =f[ f d 2 z k ff F(^_ 1 , z 2j ) det [{z^\, z^ 1 }] . (3.10) 

J.J. / J. J. ■> J Ka<2N;l<b<N L J 

fc=i J c j=i 

Here the last determinant is simply the Vandermonde determinant, and thus we have 
arrived again at eq. ( 13. ip . 

All of the above arguments can be generalised, by including an arbitrary observable 
that is symmetric under the exchange of all the eigenvalues. Perhaps the simplest 
example, which we will also encounter in the next section, is a factorising operator, 

2N 

f(x,z,z*) = nx)Hf(z 3 ). (3.11) 

3=1 

The individual factors do not affect the symmetry arguments above and we obtain 

1 » 2N r N 

(f(x,z,z*)) 2N+x -- — — I dyX(h(y)f(y))xl[ I d 2 z k f(z k )Y[F(z 2j _ ljZ2j ) 

3=1 



Z 2N+X N\ 



fc=i 



xA x+2N (y,{z}) (3.12) 

for general expectation values. An explicit example for such an operator is the 
characteristic polynomial. This result can be generalised to non-factorising observables 
symmetric in the eigenvalues. Since the monomials in the traces of a matrix can be traced 
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back to products over characteristic polynomials eq. (I3.12p is true for all symmetric 
polynomials in the eigenvalues. With WeierstraB' approximation theorem all symmetric 
functions are built of polynomials in the traces, and as a limiting case one can also 
construct distributions like the Dirac distribution. This means that all observables 
symmetric in the eigenvalues fulfil a similar equation as ( 13. X2I) in a weak sense. 

As a further remark the symplectic non-Hermitian ensembles with /3 = 4 0, [6] are 
already of the factorised form eq. (13. ip from the onset, with 

F^= A \z u z 2 ) = i(z 2 - Zl )w(z u z 2 )5 2 (z 2 - z{) (3.13) 

and we always have \ = due to symmetry. The symmetric real weight w(zi, z 2 ) can be 
found in [6] and [9] for the Gaussian Ginibre and chiral classes respectively. Therefore 
all statements we derive from the form of eq. (13.11) automatically hold true for these 
symmetry classes as well. The factorisation thus unifies the non-Hermitian ensembles 
for P = 1 and 4; in fact, this was already pointed out in [HI [21]. In [H], this was found 
in a different way without using factorisation. 

4. Integral representation of skew-orthogonal polynomials and their 
Cauchy transforms 

In this section we will derive integral representations for the SOP for general weight 
functions, using the results from the previous section. For this we will only need the 
result for an even total number of eigenvalues (i.e. x — 0)- 

All the matrix or complex eigenvalue models introduced previously can be solved 
for all eigenvalue density correlation functions in terms of the following skew-symmetric 
kernel 

2N-1 



where 



lC 2N (z 1 ,z 2 )= 2_^ A k l Pk(zi) pi(z 2 ) (4.1) 

k,l=0 

A ki = 2 d 2 z 1 d 2 z 2 F(z 1 ,z 2 )p k (z 1 )p l (z 2 ) . (4.2) 



In fact, the kernel is only a property of the measure F(z 2 j-±, z 2 j) and not of the particular 
choice of the polynomials {p k (z)}; in [27] these were chosen to be monic. For an odd total 
number of eigenvalues a similar representation to eq. (14. ip holds, but with a modification 
to the last row and column of the matrix A; see [28] and [30]. Here we will choose the 
polynomials Pk{z) to be skew-orthogonal with respect to the following anti-symmetric 
scalar product 

f(zi) g(zi) 



f(z 2 ) g(z 2 ) 



(4.3) 



(f\g) = -(g\f) = j d 2 Zl d 2 z 2 F( Zl , z 2 ) det 
c 2 

defined for two functions f(z) and g(z) that are integrable with respect to the weight 
functions contained in F(z\, z 2 ). This includes the particular function g(zi, z 2 ) from eq. 

(EZD. 
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Our skew-orthogonal polynomials qk(z) are defined to satisfy 

\Q2k\Q2l+l) — hkSkl , 

(q2k\q2i) =0 = (fcfc+il&H-i) VA;,/>0 , (4.4) 

where the hk > are their positive (squared skew) norms, see eq. (14.161) . This leads to 
a block diagonal matrix A = diag (hoe, . . . , h^-\e) that can be easily inverted, where e 
is the anti-symmetric 2x2 matrix with elements e\ 2 = 1 = —e 2 \. The kernel can be 
written as a single sum in terms of the SOP: 

N-x 1 

K-2n(Z\, Z 2 ) = Y^ 7- (g2fe+l(2l)g2fc(^2) - q2k+l(z2)q2k(Zl)). (4.5) 

^^ fife 

The kernel for 2iV + 1 contains the same SOP plus a correction term, see also [30] for 
the Ginibre ensemble and [1] for the GOE. 

4-1. Skew- orthogonal polynomials 

After all this preparation we come to our second result, an explicit representation of the 
SOP. They are given in terms of the following expectation values 

g 2n (*) = (det(z-J)) =(f[(z-Zj)\ (4.6) 

W / 2n 

for the even polynomials, and 

q 2n +i(z) = (det(z-J)[TiJ + z + c\) = (det(z - J)Trj\ + (z + c) q 2n (z) 

I 2n \ I 2n 




in 



2n 

Zi + Z + C 
1 



(4-7) 

2d 

for the odd polynomials, which are both expectation values over an even number of 
eigenvalues In, n > 1. For n = we simply have q , o(^) — 1 an d qi(z) = z + c, by 
definition. It is easy to see by taking large arguments that these representations are in 
monic normalisation, viz q n (z) = z n + 0(z n ~ 1 ). Similar expressions hold in terms of the 
matrix T> from eq. (12.41) for the chiral model (see subsection 15.21 for more details), whilst 
the representations given in terms of squared eigenvalues A| = Zj are identical. Eqs. 
(14. 6 p and (14. 7p were also shown for particular non-Hermitian ensembles in refs. [T2"l[T3"] . 

The set of odd polynomials is not unique because of the anti-symmetry of the skew 
product (14. 3p . as the arbitrary constant c times the even polynomial drops out. In most 
of the following we will keep the constant c^O though. 

The proof of the first integral representation eq. (14.61) goes as follows. We can write 
the product times the Vandermonde determinant of dimension 2n as a Vandermonde 
determinant of dimension 2n + 1, and so 

1 2n „ n 

q2n(z) = — I '/ d 2 z k \\F(z 2] . 1 ,z 23 ) det [{z a b ~ l }\z a - 1 ] . (4.8) 

Zo„ - LJ - / - LJ - Ka<2n+l;l<b<2n L J 

2U k=l -L 7=1 - - , - - 
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We can now apply a slight modification of the generalisation of the de Bruijn integral 
formula proved in Appendix C.2 of [29] . 

2n „ n 

{fa(Z2j-l),g a (Z2j)}\ a 



Yl d 2 z k Y[F(z 2j ^,z 2j ) det 

fc=l /, j = l - - - - 



\m(m-l)/2 n j pf 

l<a,b<2n+m; l<i<m 



f d 2 ud 2 v F(u,v)[f a (u)g b (v) - f b (u)g a (v)} \ a 





-a 



<t 



(4.9) 



it 



The overall sign can be seen by choosing [a ai ] = [ mj 2 n l m . Y ■ Here a is a constant 
matrix, which in our case in eq. (14.81) is a simple vector with m — 1. In contrast to the 
usual de Bruijn formula we integrate over In variables here instead of n, as is shown to 
hold in Appendix A| (see also Appendix C.2 in |29j). 

Denoting the basis functions of monic powers by e a (z) = z a and using the fact that 
we have equal functions f a (z) = g a ( z ) — e a -i( z ) above we arrive at 



q2n{z) 



nl 



Pf 



Z,2 n l<a,b<2n+l 



{(e a _i|e 6 _i)} e a _i(^) 
-e 6 _i(z) 



(4.10) 



It is easy to see that using the definition of the skew product in eq. (14. 3 p and performing 
one more integration we have 

(<?2n|e c > = Vc = 0,...2n, (4.11) 

because the corresponding Pfaffian vanishes. Using the linearity of the skew product 
we can deduce that the even polynomials q2n( z ) in eq. (14 .6p are skew-orthogonal to all 
polynomials of lower and equal degree. 

To prove the second integral representation eq. (14. 7p we need a further identity for 
manipulating Vandermonde determinants, 

1 ... 1 

N Z\ . . . Zn 



^2z a A N ({z}) = det 



a=l 



N-2 



,N 



JV-2 

'N 



Aiv(W) 



(4.12) 



which is proved in Appendix B We can now proceed as in eq. (14. 8p . by first 
incorporating the product in eq. (14. 7p into a larger Vandermonde determinant, and 
then applying the identity (I4.12p for 2n+l: 

.. 2n „ n 

q2n+i(z) = ^-\[ d 2 z k Y\F(z 2j -i,Z2 j ) det 

ZsOr, I Ka,o<2( 

ln k=\ J c j=l ~ ~ 

For simplicity we have set c = here as it does not affect the proof. Again applying the 
integral formula eq. (14. 9p . with a slightly modified range of indices compared with the 



{z^ 1 } z a ~ x 



„2n+l 



Z 



2n+l 



(4.13) 
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even polynomial case, we obtain 

{(e a _i|e 6 _i)} (e a _i|e 2n+ i) e a _i(z) 



Til 

Q2n+l{z) = ~= L - Pf 

Z 2n l<a,b<2n+l 



(4.14) 



(e 2n +i\eb-i) e 2n +i(z) 

-e 6 _i(z) -e 2n+1 (z) 

On taking the skew product f!4.3[) of this result it obviously follows that 

(g 2 n+i|e c ) = Vc = 0,...2n-1 , (4.15) 

this being the skew-orthogonality of the odd polynomials q 2n +i(z) given by eq. f)4.7|) to 
all polynomials of degree less than or equal to 2n — 1 . 

As a last step we will verify the coefficient of the only non-vanishing skew product 
which due to linearity and eq. ( 14. lip equals (q 2n \e 2n+ i) = (q 2n \(l2n+i) ■ To do so we 
will first determine the partition function in terms of the norms. It follows along the 
lines of eq. (14. 8p . Inside the Vandermonde determinant there we could choose any set 
of polynomials in monic normalisation, after applying the invariance properties of the 
determinant. We thus have 

2n „ n 

Z 2n =\ d 2 z k ]~\F(z 2j _ 1 ,z 2j ) det [q a -i(z b )} = n\ Pf [(g a -i|<?6-i)] 

-*--■-/ ■*■■*■ Ka,b<2n Ka,b<2n 

k=l J c 3=1 ~ ~ ~ ~ 

n-1 

= n\ 

a=0 

In the second step we applied once more the integral identity (14. 9p . with m = and the 
matrix a absent. Due to the skew-orthogonality, the matrix inside the Pfaffian becomes 
block diagonal, with the norms h k times e down the diagonal, which finally leads to the 
product of the norms. 

Using eq. (I4.10p after replacing the monic powers with polynomials q k we have 



n^. (4.i6) 



77/! 

(<?2n|<?2n+l) = ~=~ Pf 

Z 2n l<a,b<2n+l 



{(q a -i\qb-i)} (q a -i\q-2n+i) 

-(Q2n+l\Qb-l) 



= ^-l\h a = K, (4.17) 

and thus the consistency of the normalisation of our integral representations ( 14. 6ft and 
( 14. 7p with respect to eq. ( 14. 4p . This concludes our proof of the integral representations 
of the SOP satisfying eq. ( I4.4p . An entirely different derivation of the same results can 
be made by a mapping to the (3 = 4 symplectic case. When mapping our F(zi,z 2 ) 
as in eq. (I3.13P we could in principle copy the orthogonality proof from [6] where the 
representations eq. (14. 6p and (14.71) were derived for (3 = 4 . 

It is worth mentioning that the same representation for SOP holds for Hermitian 
RMT at (3 = 1 and 4 with real eigenvalues as was shown earlier in [T2J [13]. However, for 
all four cases - two Hermitian and two non-Hermitian - the jpdf and corresponding skew 
products are different. In contrast for = 2 all OP are obtained from a single relation 
as in eq. ( 14. 6p . in both Hermitian and non-Hermitian RMT [22]. The fact that the same 
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integral representation for SOP holds both in non-chiral [I2J d3] and chiral ensembles 
is straightforward in the Hermitian case. However, for non-Hermitian ensembles this 
becomes nontrivial comparing [6] vs. [9] for /3 = 4, and [TD] for (3 = 1. This is due to the 
two-matrix model structure of the chiral ensembles, where the change to an eigenvalue 
basis requires detailed calculations. 

Let us finish this subsection with some remarks. For an even number of eigenvalues 
X = the anti-symmetric kernel eq. (14.11) (or (14.51) ) can itself be expressed as the 
expectation value of two characteristic polynomials for (3 = 1 |10) ||l (and f3 = 4 [21] ) 



det(A - J) det( 7 -J)) = hN ^™^^ with A ^ 7 , (4.18) 

/ 27V A — 7 

and similarly for the chiral ensemble. This equation is valid for arbitrary weight 
functions. In fact we will partly use this relation to determine the set of odd polynomials 
eq. ( 14 .7p in section [5] below. So why are eqs. ( 14. 6 p and ( 14. 7p interesting if the kernel 
itself can be independently determined as a building block? It is because the integral 
representations we just derived, and the explicit determination of the SOP in some 
examples in the next section, complete the list of classical polynomials in the complex 
plane for the three elliptic Ginibre ensembles and their chiral extensions. 

The determination of the SOP through an ansatz, and subsequently the direct 
verification of the relations ( 14.4p for skew-orthogonal Hermite polynomials, was already 
a formidable task for the elliptic real Ginibre ensemble as can be seen from [TJ. Because 
of the non-Gaussian form of the chiral weight eq. ( 12 .9p this is even more so true for skew- 
orthogonal Laguerre polynomials. The integral representations derived here provide an 
alternative, constructive approach, leading to a new result for skew-orthogonal Laguerre 
polynomials. 

4-2. Cauchy transforms 

We now come to the definition and integral representation of the Cauchy transforms 
tk(z). It is very natural to define them with respect to the scalar product eq. ( I4.3P as 
follows: 

<ln{zi) 



t n {n) = / d 2 zi d 2 z 2 F(zi, z 2 ) det 



K — Z\ 

q n (z 2 ) 



K- Z 2 



(I,, 



1 



K — Z 



(4.19) 



C 2 

We will now show that the following integral representations hold: 

*■<«> = "- (stobjL = "" ( 5 (^y )_ (4 - 20) 

for the Cauchy transforms of the even polynomials, and 

, , , /TtJ-(k + c)\ , / TrJ \ . , . , 

'-'<"> - h "{ detfc-7) / 2 „ +2 = ft "\d^7)/ 2 „ +2 - < K + ^«M 

|| Note that the overall constant in front of the kernel has been chosen here to be consistent with the 
standard choice in eq. (j4.5[) . 
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2n+2 

Y, Za ~ (« + C) 
o=l 

2n+2 

no*- *i) 



(4.2i; 



2n+2 



for the odd polynomials. Note that the averages for £2™ (ft) and £2n+i(ft) run over 2n + 2 
variables, instead of 2n as for the polynomials q 2n (z) and (fen+i^)- This implies in 
particular that io(ft) 7^ constant, see also eq. (14.221) below. 

The correct overall prefactors can also easily be seen. From expanding the geometric 
series in the definition eq. ( 14.191) for large arguments, and using eqs. ( 14.111) and ( 14.151) 
as well as the anti-symmetry of the first non- vanishing skew product, it follows that the 
Cauchy transforms are indeed Laurent series with the following coefficients 

K „ / 1 



Wft) 



K 



^2n+l(ft) 



2n+2 
K 

-271+1 



o 



+ o 



-2n+3 



-2n+2 



(4.22) 

The form given in eqs. ( I4.20p and (j4.2ip is completely analogous to eqs. ( 14. 6 p and 
( 14. 7p . as well as to the corresponding result for (3 = 2. Let us also remark that such a 
representation was not known before in the non-Hermitian (3 = 4 symmetry class, and 
that both translate into new representations for = 1,4 in the Hermitian limit. 

We begin by proving the representation for the Cauchy transforms of the even 
polynomials. Inserting the result eq. ( 14. 6 j) into the definition ( 14.19(1 we have 

(det(2i - J)) 2n (det(^ 2 - J))ln 



t 2 n\ K ) 



d z x d z 2 F(zi,z 2 ) 



K- Z 2 



Zl 



— — / d z 2n+ id z 2n+2 F{z 2n+ i,z 2n+2 ) y j d ZkyyF{z 2 j_i,z 2 j) A 2n ({z}) 

1 — 1 J 1—1 



Ilj=l( z 2n+1 - Zj) Yl j=1 (z 2n+2 - Zj] 



~ Z 2n+2 
2n+2 



K 



Z 2n +l 



n + l)Z. 



2 n 



2n+ 2 „ n+1 

JJ / d 2 z k Y[F(z2j-i,z 2 . 

k=l I j=l 



n+1 

xE 

fc=i 

1 



A 2 n+l({z}«^2fc) A 2 n+l({z}l^ 2 k~l] 



ft- z 2k 

2n+2 



K 



Z 2 h-\ 



[n 



\)Z. 



'in 



fc=i 



n+1 



d*z k \\F\ 



Z 2 j-1, Z 2 j 



2n+2 



[n 



1)Z 



2n 



fc=l 



n+1 



det 

l<a<2n+2;l<fe<2n+l 



k- 1 } 



Za 



I d*z k ][F(z 2j _ u z 2j ) ] 



2n+2, 



3=1 



1=1 



Zl 



det(«; — J) 



(4.23) 



2n+2 
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In the first step we have simply written out the expectation value and renamed the 
additional two integration variables. The products in the numerator can be incorporated 
into a larger Vandermonde determinant. Next we can symmetrise the integrand with 
respect to an exchange of any pair of variables z 2 j, z 2 j+\ leading to a prefactor l/(n + 1). 
The resulting expression can be seen to be the expansion of a Vandermonde determinant 
plus an extra column. In the last step we use an identity that was proved in [29], see 
eqs. (3.3) vs (3.7) therdJJ deriving different representations for Berezinians 



det 

Ktt<2n; K6<2u-1 



K" 1 } 



n 



2ra i 
a>b\ 



z b) 



n£i(«-*i) 



(4.24) 



This can be used to express the Cauchy transform as an expectation value, after 
providing the correct normalisation factor from eq. (14 . 1 6 j) in the last step. 

The proof for the odd Cauchy transforms follows along the same lines. For simplicity 
we set c = here, which can easily be reinstated at the end: 



^2n+l(^) 



1 
^2n 



2ra+2 

n 



fc=i 



c 

2d 



n+1 

d 2 z k Y[F(z2j-i,z 2j ) 



^2/l 



Ezn . 

i=\ z i ~r z 2n +i 



\[j=l{ Z 2n+l 



K — z 2n+2 



Z)j=l Z l + Z 2n+2 ) ri 7 =l(^2 



2n+2 



Zj) 



z 2n+l 



K 



A 2 n({4) 



1 



(n + l)Z 



2n+2 „ n+1 

-] d 2 z k \\F{z 2j 

2n k=i i j=i 



-li z 2j) 



n+1 

><E 

fc=i 
1 



^2n+l({ z }l^2k) ^2n+l{{ z }l^2k-l, 



K - Z 2k 
2n+2 



K 



z 2k-l 



[n + l)Z 2n 



ln+2 „ n+1 

nd 2 z k ]~\ F(z 2 j-x,z 2 j) det 

/ J-- 1 - J J l<a<2n+2;1 

fe=i i j=i 



<fe<2n 



l^" 1 } 



2n+l 



K 



z a 



2n+2 



n + l)Z. 



2n 



fc=l 



n+1 



d 2 z k \\F{z 2j _ 1 ,z 2j ] 



E 



2n+2 
=1 



Zl 



^jnSTte-^) 



h„ 



TrJ 



c 



i=i 



n^(« - zi) 



(4.25) 



Here we included the product into the Vandermonde determinant as before, as well 
as the additional sum leading to the modified Vandermonde determinant A defined in 
the right-hand side of eq. (j4.12p . In the last step we simply need a slightly modified 



version of the identity eq. (14.241) which is derived in Appendix C The inclusion of the 
arbitrary constant c 7^ follows simply by shifting k — > K + c in the numerator but not 
in the denominator. This concludes the derivation of all the integral representations of 
the SOP and their Cauchy transforms. In principle the simple expectation values in 

If Note that we order products here so that there is no sign in eq. (|2.8jl for the Vandermonde 
determinant. 
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eqs. ( 14. 6 p and (J4.7P as well as eqs. (j4.20p and (j4.2ip could be computed explicitly using 
supersymmetric vectors depending on ordinary variables and Grassmannians. In the 
explicit examples given in the next section we shall give the resulting SOP only. 

5. Examples for skew-orthogonal polynomials 

In this section we will give three examples of skew-orthogonal polynomials in the complex 
plane: Hermite, Laguerre and Laguerre with mass terms. Although the first of these 
were already known, our derivation is new. The second and third are new examples. 

In principle there are two different methods. In the first of these we directly use the 
integral representations; for the even polynomials these are the expectations of a single 
determinant, eq. (14. 6p . and for the odd polynomials the expectations of a determinant 
multiplied by a trace, eq. (14. 7p . Both can be calculated in one step by computing the 
expectation of the product of two determinants (which is proportional to the kernel) 
and either taking limits or differentiating, and using the fact that the determinant is the 
generating functional of all independent matrix invariants. The expectations can then 
be computed using Grassmannians, and because this was already explicitly done in [10] 
we can be very brief here. 

The second method follows the general setup outlined at the start of section HJ 
Given the kernel in terms of general polynomials eq. (14. ip . and choosing them to be 
skew-orthogonal eq. (14. 3p . the individual polynomials can be "read off" from the kernel 
in eq. (I4.5P by differentiation (or taking limits): 

q2n{Z) = K (271 + 1)! da*** /C2 " +2( ^ ^ = K uToo M 2n+1 

1 d 2n 
q 2 n+i(z) = ~h n J C 2 n+2(u,z) + c q 2n (z) . (5.1) 

(2n)\ ou zn u=o 

This is possible whenever the kernel has already been independently determined, e.g. 
by the above procedure detailed in [10] (see also [27] for another method). In addition 
the norms hk can be read off from the kernel as the leading coefficients. 

Of course both methods lead to the same answer. In the third example the kernel 
including the masses as well as the partition function itself have not previously been 
computed explicitly and so this also constitutes a new result. 

5.1. Example I: skew- orthogonal Hermite polynomials 

In general the calculation of the expectation of a single determinant (or the product 
of two determinants) is straightforward, even without switching to an eigenvalue basis: 
we express the determinant as an integral over anti-commuting (Grassmann) variables, 
and then the Gaussian random matrices can be integrated out. After a Hubbard- 
Stratonovich transformation the anti-commuting variables can also be integrated out. 
Because the procedure was carried out and explained in detail for two determinants 



(5.2) 



toM = to -! = ^ = ( i ) #»(-£=). (5.3) 
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with Nf = in this model in [TD] we only quote here the result for our first example, 
the expectation with respect to model eq. (I2.3ll + I: 

n i / \ / 

det(z - J) det( U - J)) = JVI J>< £ m H k (-j=) H„ f-|= 

where r is the non-Hermiticity parameter, and the Hk(z) are the standard Hermite 
polynomials. Hence, for the even polynomials we can simply project out the second 
determinant to give 

det(z — J) det(w — J)) /r\ k 

«->oo u 2k \2/ ZK \\/2t 

Here we used the following result to calculate the single term in the double sum that 
survives the limiting process: 

u-kx) u \a/ \ a / 

This equation also implies that the even polynomials eq. (I5.3P are in monic normalisation 
as they should be, starting with qo(z) = 1. Alternatively we could of course have 
differentiated eq. ( 15.2ft N times with respect to u. 

For the odd polynomials, we use the fact that the determinant is the generating 
functional for symmetric functions, and in particular for the trace: 

i d N - x 

— r— r detfuljv - J) = -TrJ, (5.5) 



(JV-1)! du 1 



u=0 



where J is an N x N matrix, and N > 1. Applying this to eq. ( 15. 2 p on the left-hand 
side allows us to obtain q 2 k+i(z) from eq. (14. 7p . for k > 1: 



1 & 



2k-l 



m+l{z) = ~ j2J—Ty. d^\ det{z ~ J) det{u ~ J) / 2fc 



+ {z + c)q 2k {z) 

u=0 



(2*) ( «/I ) (r + !)#»_! f -jjL ) + (z + c) ( I) # 2fc 



2; v ' \V2tJ v 2/ VV2T 

k+l / \ fr-i / \ k 

T 2 „ I Z \ , T 2 I Z \ T 



-H 2k+l [-= )+ 2k — -H^ -= + c-jrflak -= . (5.6) 



2 fc+ 5 Vv2r/ 2 fc 2 \V2tJ 2 Vv2r, 

In the first step only two terms survive the differentiation after setting u = 0; we used 
the following properties of Hermite polynomials in addition to eq. ( 15. 4ft 

— tf n (.) = 2^U (5.7) 

if n+ i(z) = 2zH n (z) - 2nH n _ x {z) , for n > 1 , (5.8) 

as well as the recurrence relation to simplify eq. (I5.6P in the last line. This form makes 
it more transparent that the arbitrary addition of cq 2 k{z) is the only even Hermite 

+ The double sum can be simplified by using the Christoffcl-Darboux identity [TU] . 
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polynomial appearing in this example. From eq. ( 15. 4p it also follows that q 2 k+i(z) is in 
monic normalisation, and for k = we have qi(z) = z + c by definition. Defining 

a M = (£)*M^) (5 ' 9) 

we reobtain the following final simple result from [7] 

q2k{z) = C 2k (z), 

q2k+i{z) = C 2k +i(z) - 2kC 2k -i(z) - cC 2k (z) . (5.10) 

The norms hk = 2(r + l)v27T (2fc)! can be determined either by direct calculation 
of the scalar product eq. (14.41) of the SOP which we just obtained, as was done in [7J, 
or by computing the partition functional 

Let us emphasise that in our derivation the skew-orthogonality of the polynomials 
is automatically satisfied due to their general integral representation eqs. ( 14. 6ft and ( 14. 7ft . 
This is in contrast to [7J, where the skew-orthogonality was explicitly verified for the 
weights eq. ( 12. lip in the complex plane. 

5.2. Example II: skew- orthogonal Laguerre polynomials 

In this subsection we turn to entirely new expressions for skew-orthogonal Laguerre 
polynomials for our chiral model eq. (12. ip . We will start with the so-called quenched 
case (Nf = 0) and then add mass terms in the next subsection. 

Expressed in terms of the 2n eigenvalues Zj of AB T these are the same as before; 
however, the matrix T> in eq. ( 12.41) also has v generic zero eigenvalues, and so we repeat 
the integral representations eqs. (14.61) and (14.71) for completeness, and also to make 
contact with [10] . For the even polynomials, we have 



q2n{z) = ^ E ( det(v^l in+v-V)) =(det(zt 2n -AB T )\ (5.11) 

Z V I A \ / An+v \ I In 



and for the odd polynomials 

q2n+l(z) = -^ (det(y/zt 4n+u - V) 



-Tr£> 2 + z + c 
2 



ln.+v 



det(zt 2n -AB T )TiAB T ) + (z + c) q 2n (z) . (5.12) 

/ 2n 

The starting point for what we need for our calculations, namely the expectation of 
two determinants, was again given in detail in [10] and thus we merely state the result: 

/ detail 2n - AB T ) det(wl 2 „ - AB T )\ 

\ I 2n 

- ( 2n)! (2„ + vW M« t (I) " t (S^f ^ (^) H (^) , (5,3) 

where we recall the notation ( 12. 2p . 

* The lower order terms from combining eqs. (|4.5j) and (|4.18|) will provide ratios of norms Hn /hk and 
thus the /c-dependence only. 
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We thus obtain the even polynomials by simply projecting out the second 
determinant 

q 2k (z) = lim — : (det(zl 2k ~ AB T ) det(ul 2fc - AB T 

= (^Vf(2fc)!^(^-). (5.14) 

Here we have used the following relation for the Laguerre polynomials 

1 ra\ (-1) N 

lim — L n{-) = T7TF' ( 5 - 15 ) 

u->co u" \aJ TV! a" 

which also confirms that the even polynomial is properly normalised. 

For the odd polynomials we again need to take derivatives as in eq. (15. 5 1) to obtain 

for k > 1 

^^( det(zt 2k - AB T ) det{ut 2k - AB T )) 

u=0 



««*+lW = - (^3T)T Ih^i^ 12 * " ABT) det(Ml2fc " ABT) /2k 

+ (z + c)q 2k (z) 

-.2 



(V,_)-'(2*)! (2k + „) ( 2 ^ t (^-i-) + (l + |)L Vl (^)) 



+ (z + c)(V,_)* (2*)! L^ 



Z 



4yU 2 77_ 

(4A_) 2fc+1 (2* + 1)! L 2 Vi {t^A + c '(V^) 2fc (2fc)! ^ 



J2k \4^ 2 rj. 

+ (2fc + ^A+ftVV) 5 *- 1 (2A;)!L 2 V 1 (^") > ( 5 - 16 ) 

where the new arbitrary constant 

c' = c+(4/iV)(4fc 2 + 4A; + l + (2A; + l)i/) (5.17) 

now depends on k, v and /x. The above result was obtained after using the corresponding 
relations for Laguerre polynomials: 

d n ~ x 
J-^I L n(*) =(-l) n (z-(n + u)) (5.18) 

(n + l)L u n+1 (z) = {2n + v + l-z)L v n {z)-{n + v)L u n _ 1 {z) ,for n>l. (5.19) 

It is easy to see that the polynomials are again monic, due to eq. (j5.15p . This once more 
fixes q\ (z) = z + c. We can now define 

Cl{z) ee (4/iV) fc k\ LI (AA (5.20) 

allowing us to write 

q2k{z) = +C% k {z), 

fcw-i(*) = - C» k+1 {z) + (1 + /i 2 ) 2 (2fc)(2fc + v)C^_ x {z) + JC» k (z) , (5.21) 

giving the new skew-orthogonal Laguerre polynomials up to an arbitrary constant. The 
final result compares with the similar form of eq. (j5.10p . 
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For the norms we find h k = 8n(4fi 2 )(2k)\ (2k + v)\ (Afi 2 n + ) 4k+u+1 where the k- 
dependence again follows from the ratio of the norms h N /h k , see eqs. (14.181) and (14. 5p . 
whereas the overall constant factor can be deduced from the partition function, see eq. 
(3.46) in [IT] , and taking the ratio for consecutive values of N. 

5. 3. Example III: inclusion of mass terms in the chiral model 

Our third example gives the SOP again for weights including Nf > mass terms, which 
is also called the unquenched case. We will exemplify this using the chiral model eq. 
(12. ip where such terms are more common due to applications in QCD. However, the 
same insertion of mass terms can be done in the non-chiral model eq. (12.31) following the 
same lines. 

Our main point here will be to express the SOP for Nf > in terms of the SOP for 
Nf = (and the corresponding kernel), which we have already calculated. To indicate 
which polynomials we are referring to we will use a superscript, as in q k f (z), and 
correspondingly for the kernel and expectations. 

Our derivation relies heavily on [21] where all the expectation values of products 
and ratios of characteristic polynomials (or determinants) have been expressed in terms 
of Pfaffian expressions of matrices containing a small number of building blocks; in our 
case these building blocks will be the quenched (Nf = 0) SOP and kernel. 

To begin we first express the unquenched integral representations eqs. ( 14. 6ft and 
(14. 7p in terms of ratios of quenched expectations. For the even polynomials we have 

<lin f \z) = ( det(zt 2n - AB T )) {Nf) (5.22) 



2v 



( det(zt 2n - AB T ) \f f L x det(m}t 2n - AB 1 



(o) 



2/i 



(0) 



U%det(m}l 2n ~AB T )) 

1 J I In 



and similarly for the odd polynomials we have 

( det(zl 2n - AB T ) TiAB T U% det(m}t 2n - AB T )) (0) 
&\(z) = ± -^ ^l + ( Z + c)&'\z). 

{U%det(mjt 2n -ABT)j 2n 

(5.23) 

We will now give all the building blocks for these expressions. The first building block 
in the denominator, the expectation value of the mass term, simultaneously provides us 
with the massive partition function itself: 

7 {N ? ) (Sm\) Nf I Nf x (0) 

""if !) = IHf ( II det Kl- - AB-))[l . (5.24) 

^ch2N /=1 /=1 

The masses to the power of z/, the number of generic zero eigenvalues, of course cancel 
in the ratios for the q k i (z) above. Using the results from [21] and expressing the 
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expectation values there in terms of our quenched kernel and even SOP we obtain: 
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X Pf 
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where we have to distinguish even and odd numbers of flavours Nf 
the norms in the prefactor is equal to unity when the upper limit is N 
to [21] we have used the following identity 

?2m(™« 



N f odd (5.25) 

The product over 
1. Compared 
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(5.26) 



which can be easily seen by adding multiples of the first row and column to the remaining 
rows and columns, in order to eliminate the leading SOP in the kernels and hence shifting 
their index down by two. This result for the partition function (or expectation values of 
characteristic polynomials) precisely equals the corresponding result for (3 = 4 in [24p. 
The even polynomials now easily follow from eq. ( 15.251) . by choosing one of the 
masses to be the argument. We obtain 
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(5.27) 



for Nf even. Here and in the following we suppress the indices of the Pfaffian which run 
from 1 to Nt in both the even and odd cases. For Nt odd we obtain 
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(5.2* 



Next we determine the massive kernel using eq. (14.181) 
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Notice a typo in 



U%det(m 2 t 2n -ABT)) {0) 

3 J I 2N-2 

in eq. (2.8) compared to the correct Theorem 1 in eq. (3.1) there. 
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(5.30) 



Here the mass dependent inverse norm l/h N £ l has been eliminated using the following 
identity, leading to a shift in the index of the kernels in the denominator by +2. 
Following eq. (I4.16P we can write 
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N f odd. 



Likewise for Nf odd we have 
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(5.32) 



Let us pause with a few remarks. Following [TT] these expressions for the massive 
kernel determine all massive eigenvalue correlation functions for 2iV, in terms of the 
known quenched kernel and the quenched even SOP that were given in Example II 
above. In particular it is transparent that even for finite N the unquenched kernel 
(when properly normalised by the massive weight) is given by the quenched kernel plus 
some correction terms. The same structure thus prevails for all unquenched eigenvalue 
correlation functions. 

For 2N + 1, correction terms to the massive kernel will also include the massive 
SOP q 2N i (z), when following e.g. [30]. Thus all massive eigenvalue correlation functions 
follow in this case as well. 
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As the final step we will give the massive odd SOP q 2N f +1 (z). Here we will follow 
eq. (15. ip and determine them from the kernel, rather than eq. (15. 5p . As an aside, above 
we could have alternatively determined the kernel first and then the even SOP from eq. 
(15. ip as well. A slight generalisation of eq. (15. ip reads 
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u=0 



where k > 0, and the a; are some arbitrary constants. It is easy to see that we only 
get a non-vanishing result when 2n or 2n + 1 of the derivatives act on the kernel and 
not the prefactor. This is true because the function in the brackets is a polynomial of 
order 2n + k + 1 in the variable u. Hence, the differentiation yields the coefficients of 
the monomials of order 2n and 2n + 1 in the variable u of the kernel K, 2n f +2 . 

In choosing k = Nf and the ai = m 2 we can use this relation to cancel the 
factor rif=i( u ~ m t) m the denominator of eq. ( I5.30p that would otherwise have to 



be differentiated as well. We thus obtain from eq. (I5.33P and (15.301) that 
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K-2N+N f { rn }i m g) 

for Nf even. 

Here we have pulled the derivatives inside the Pfaffian leading to the quenched SOP of 
shifted odd index, and used eq. (15.311) and the identity corresponding to eq. (I5.26P for 
the odd polynomials. For odd Nf we obtain 
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+ c'q ( 2N f, ( y z) for N f odd. 



(5.35) 



The constants c in the Pfaffian in the numerator can be absorbed into the even 
polynomial c'q 2N f (z), which can be seen as follows. Just as the determinants of two 
matrices that only differ by a single row (or column) can be added, a similar statement 
holds for Pfaffians: due to linearity the Pfaffians of two anti-symmetric matrices that 
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only differ by a single row and its transposed column can be added. We can thus split 
off the c-part from the Pfaffian above to obtain 
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K-2N+N f +1 ( m /' Z ) 

which is proportional to the numerator of the even polynomials with odd Nf in eq. 
(15.281) . Thus the final result for the odd polynomial with odd Nf is eq. ( 15.351) with c = 
and d —> c . This ends our third example for the SOP and kernel including masses. 

Similar expressions could be given for the non-chiral model eq. (I2.3p . as well as for 
the Cauchy transforms of the unquenched SOP. 



6. Conclusions 



In this paper we have completed the analysis of the set of (skew-) orthogonal polynomials 
in the complex plane that apply to the chiral extensions of the three elliptic Ginibre 
ensembles. By constructing an explicit integral representation we found a new set of 
skew-orthogonal Laguerre polynomials in the complex plane which provide an alternative 
method of solving the chiral ensemble with real asymmetric elements. Our integral 
representation is also valid for the real elliptic Ginibre ensemble; in fact, it holds for 
arbitrary weight functions g and h in these two classes. Furthermore we also gave a 
new integral representation of the Cauchy transforms of these polynomials which holds 
not only for the two symmetry classes with real matrix elements (/3 = 1) but also for 
quaternion real matrix elements (/3 = 4). 

An important ingredient for our results was a proof that the probability distribution 
in the partition function factorises for /3 = 1. This offers another unifying view of the 
non-Hermitian /3 = 1 and f3 = 4 symmetry classes, both chiral and non-chiral. 

There are many more non-Hermitian ensembles, in addition to the three Ginibre 
classes and their chiral counterparts, all six of which have now been solved. It thus 
remains an open question whether corresponding sets of orthogonal or skew-orthogonal 
polynomials exist for the other ensembles. It is possible that, just as in the real case, the 
known polynomials also apply to some of these other non-Hermitian symmetry classes, 
once a complex eigenvalue representation has been found for them. 

As an application of our results we have shown how to construct the skew-orthogonal 
polynomials and the kernel when including Nf characteristic polynomials or mass terms 
into our models. These constitute the building blocks for the massive partition function 
and eigenvalue density correlation functions. Consequently this will allow us to study 
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the complex Dirac operator spectrum for Quantum Chromodynamics with two colours 
and non-vanishing quark chemical potential, both in the low and high density phases. 
The study of the large- N limit needed for this partly follows from the known quenched 
Nf = results and is left for future investigations. 
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Appendix A. Generalisation of the de Bruijn integral formula 

In this appendix we slightly generalise the standard de Bruijn integral formula that 
reads 
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Here w(z) is a weight function in the complex plane and / and g are functions such that 
the integrals exist. The proof is usually done by a Laplace expansion into 2x2 blocks 
that each depend on a single variable z\. 

If we start out with 2n instead of n integrations over a product of an anti-symmetric 
weight F(u, v) and let / and g depend on different variables, we have on the left-hand 
side 
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(A.2) 



Here (— ) <J is the sign of the permutation of the 2n variables, and the sum is over all 
(2n)! permutations which satisfy the restriction <r(l) < er(2) < . . . < o~(2n). Note that 
each pair {z 2 j-i, z 2 j} only appears in one subdeterminant. This gives the Pfaffian as a 
result (see e.g. P for a definition). 
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Appendix B. Relation to a modified Vandermonde determinant 
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In this appendix we prove an identity related to Vandermonde determinants, which is 
needed to derive the integral representation eq. (J4.7P for the odd polynomials g 2 n+i(z)- 
For completeness we repeat the relation eq. ( 14.121) which is to be shown here, 
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Using the second representation from eq. (12. 8p . A N ({z}) = deti< 0] 6<jv[^a 1 ], one can see 
that the modified Vandermonde determinant in eq. (IB.lj) has a mismatch of 1 in the 
powers in the last row compared to the Vandermonde determinant. 

Our proof uses that An({z}) is simply the coefficient of power u N_1 in a Leibniz 

expansion of the Vandermonde determinant Ajv+i({z}, u) of size N + 1 with respect to 

the last column in the extra variable u. This term can be singled out by a differentiation, 
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On the other hand, using again eq. (I2.8P that A N ({z}) = Ylj>k 
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Zk), one can write 
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Combining eqs. (1B.2P and ( IB.3P we obtain the result ( IB.ip . 



Appendix C. Cauchy-type identity for the modified Vandermonde 
determinant 



In this appendix we prove the following identity, 
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which is a slight modification of identity ( 14.24[) with a mismatch by one power in the last 
but one column. In fact we will use the identity ( I4.24p to prove the above. Expanding 
the left-hand side with respect to the last but one column we have 
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[T,k=l Z k -«JIIi>i (^-^) 
= T-r2n+2/ x ' (^-2) 

IL=i («-*i) 

In the second step we used the identity (j4.24p and multiplied by unity to complete the 
product in the denominator. For the numerator we obtain the modified Vandermonde 
determinant A from eq. (j4.12p and a proper Vandermonde determinant, both of size 
2n + 2, after resumming the expansion. Writing out explicitly A2 n +2 from the left-hand 
side of eq. fIB.ip yields the right-hand side of our identity (IC.ip . 
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